q-Fermionic Numbers and Their Roles in Some Physical Problems 



R.Parthasarathy 1 
The Institute of Mathematical Sciences 
C.I.T. Campus, Tharamani Post 
Chennai, 600 113, India. 



Abstract 

The q-fermion numbers emerging from the q-fermion oscillator algebra are 
used to reproduce the q-fermionic Stirling and Bell numbers. New recurrence 
relations for the expansion coefficients in the 'anti-normal ordering' of the q- 
fermion operators are derived. The roles of the q-fermion numbers in q- 
stochastic point processes and the Bargmann space representation for q- 
fermion operators are explored. 
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q-deformed Stirling numbers were introduced by Carlitz, Gould and Milne 
[1]. Such numbers are encountered in [2] in the 'normal ordering' of q- 
deformed boson oscillator creation and annihilation operators [3]. Recently, 
Schork [4] considered generalized q-Stirling numbers and obtained useful 
properties. He further introduced unsigned q-deformed Lah numbers. Subse- 
quently, he [5] considered q-fermionic Stirling numbers which are encountered 
in the 'normal ordering' of q-fermionic oscillator creation and annihilation 
operators introduced by the author and Viswanathan [6]. Katriel [7] has de- 
rived q-Dobinski formula for q-bosonic Bell number by using q-boson coherent 
states. In contrast to this semi-classical derivation, a probabilistic derivation 
of the q-Dobinsky formula for q-bosonic Bell number has been obtained by 
the author and Sridhar [8] by considering q-stochastic point processes. 

The Stirling numbers are crucial in many combinatorial problems [9]. 
Milne [1] in his studies on generalized restricted growth functions obtained q- 
Stirling numbers (bosonic), q-Dobinski formula and q-Charlier polynomials. 
It is very intriguing that these numbers occurring in combinatorial problems, 
arise naturally in undeformed and q-deformed harmonic oscillator algebra of 
the creation and annihilation operators. 

In this paper we consider q-fermionic numbers in more detail and obtain 
q-fermionic Stirling and Bell numbers explicitly. Further we consider the 
'anti-normal ordering' of the q-deformed bosonic and fermionic operators 
and obtain new recurrence relations. Using their Fock space states, we give a 
meaning to these numbers as expressing powers of q-bosonic and q-fermionic 
numbers in terms of 'raising factorials'. This is the counterpart of the q- 
Stirling numbers of the second kind in expressing the powers of q-numbers 
in terms of 'falling factorials'. 

In the subsequent part of this paper, we attempt to give a probabilistic 
interpretation of the q-fermionic Stirling numbers of the second kind by con- 
sidering q-stochastic point processes. This introduces q-product densities, 
a generalization of the concept of product densities introduced in 1950 by 
Ramakrishnan [10] in his study of cosmic ray cascades. Then, we give a 
Bargmann space representation of the q-fermion operators as multiplication 
by and q-differentiation with respect to quasi-Grassmann variable leading to 
differential equations involving q-differentiation on spaces of entire functions 
of quasi-Grassmann variable. 
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2.q-fermionic numbers 

Macfarlane [3] and Biedenharn [3] introduced q-boson oscillator algebra 

as 

aa) - v^a f a = q~ N/2 ; [N, a] = -a , [N, a f ] = a f , (1) 
where q > 0. By making a transformation 

A = q N/i a ; A^ = a)q N '\ (2) 

one obtains 

A4 f - qA ] A = 1; g > 0, (3) 
and the associated q-bosonic number 

1 - a n 

[n\ b = (4) 

In the above q is strictly positive. The author and Viswanathan [6] , proposed 
a non-trivial q-fermion oscillator algebra as 

ff ] + y/qf*f = <T" /2 ! W /] = [N, / f ] = / f , / V , (/ f ) 2 ^ 0(5) 
where g > 0. By making a transformation 

F = q m f . F i = pq N/ \ (6) 

one obtains 

FF f + =1 ; F 2 ^, (^ f ) 2 7^ 0, q > 0, (7) 

and the associated q-fermionic number 

In], - 1+7^' (8) 

In (5) to (8), q is strictly positive. Thus, it is not correct to replace q by 
—q in the q-boson algebra to obtain q-fermion algebra (7), as q in (1) to (4) 
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is strictly positive. Nevertheless, for mathematical expressions such replace- 
ment may be carried out. In what follows, we shall denote (8) by q-fermion 
number. This definition is of fundamental importance and is quite different 
from [n\ h . In the limit q — > 1, [n] b — > n while [n]f — > \{\ — (— l) n ) taking 
values and 1 for n even and odd. The properties of the q-fermion num- 
bers associated with (5) (which are different from (8)) have been studied by 
Narayana Swamy [11]. For q < 1, the q-fermion numbers (8) never go beyond 
1 for any value of n and for n — > oo, it asymptotically approches 0.5. On the 
other hand the q-boson numbers (4), become the usual numbers when q — 1, 
and for q < 1, as n — > oo we have, [n]b asymptotically goes to 1/(1 — q). For 
q > 1, we have for q-fermion numbers, [0]/ = , [1]/ = 1 and [n]f > for n 
odd, < for n even. 

The q-fermions described by (5) or (7) with the q-numbers (8) are different 
from the k-fermions introduced by Daoud, Hassouni and Kibler [12]. The k- 
fermion algebra is a non-Hermitian realization of the q-deformed Heisenberg 
bosonic algebra with q being a root of unity and satisfy /± = and / + ^ /_' , 
except for k = 2, for which they become ordinary fermions. On the other 
hand, q-fermions admit q real or complex and only when q — 1, they become 
ordinary fermions. 

3. Normal Ordering of q-fermion operators 

We wish to evaluate (F^F) r using (7). It is straightforward to expand, 

(F f F) r = J2^ r s (F^) s F s , (9) 

s=l 

and find a recurrence relation for T r s . From (9), it follows 

(irt^r+i = Y^r s {F x ) s F s F^F. (10) 

s=l 

From (7), we have 

F s F f = [sj/F 5 - 1 + (-l) s q s F^F s . (11) 
Using (11) in (10) and noting from (9), = , = 0, we find 

r s +1 = (-lrV-^ + W/JT, (12) 
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the desired recurrence relation for T r s coefficients in (9). The recurrence 
relation (12) with T\ = 1 is the same as that of the q-fermionic Stirling 
numbers of the second kind [5]. 

The q-fermionic Bell number introduced in [5] is 

BP = ][>;. (is) 

s=l 

An attempt along the lines of [7] for q-fermionic Bell number runs into dif- 
ficulty, q-fermion coherent states have been constructed in [13] using 'quasi 
Grassmann' variables ip. It is to be noted that the replacement of q by — q 
in the q-boson coherent states will not give q-fermion coherent states. These 
two coherent states are structurally very different. As rjrij} + ipip^ = 0; ip 2 ^ 
0, (V^) 2 7^ 0, it is not possible to use the analogue of \z\ = 1 here. So, we 
take the matrix elements of (9) between q-fermion Fock space states [13] \n > 
with n > r and use F\n >= \J[n]f\n — 1 > ; F^\n >= \J[n+ l]/|n + 1 > to 
arrive at 

which can be verified explicitly using (8) and (12). Multiplying (14) by A n 
and summing n from 1 to oo and then setting A = 1, we obtain 

B</> = (e««(l))-' f (15 ) 

where 



) = E 



oo x n 



^0 W 

(15) is the q-fermionic Dobinski formula. 

Some of the q-fermionic Bell numbers are: 

& = 1, 
B { 2 f) = l-g, 
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$P = l-q-q[2} f -q\ 

B { / } = l-9-9[2]/-?[2]J-9 S -9 S [2]/-9 3 [3] / + « Z 6 , 
$P = l-g+(-g-g[2] / -g[2]J)([2] / + g 2 ) 

+ (V-9 3 [2]/-? 3 [3]/)([3] / -9 S ) + ? 8 [4]/ + ? 10 . (16) 
In the limit q — 1, we have 

= 1 ; M /} = 0, 

£(/) = (_!)»■ if r = (mod 3) 
(-l) r+1 i/r = l(mocZ3) 
if r = 2 (mod 3). (17) 

The results (17) have been obtained by Wagner [14] in his study of generating 
functions for some well known statistics on the family of partitions of a finite 
set. 

The q-fermionic Stirling numbers of the first kind are introduced by ex- 
pressing the inverse of (9), as 



(i?t)r F r = J2S r s (F^FY 



s=l 



The matrix elements of the above between q-fermion Fock space states \n > 
,(n>r) give 

- (18) 

Using [n]f — [r]/ = (— l) r q r [n — r]/, we have the recurrence relation 

Sl +1 = (-lyq-rs^-Wfi-iyq-rS:. (19) 
Similarly, the q-fermionic unsigned Lah numbers introduced as 

[r + n-l] f \ " [r^ 
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have the recurrence relation 

d +1 = (-l) n+a - 1 g n+a - 1 £y_ 1 + [s + n]fC?. (21) 

These mathematical results reveal the feature of obtaining them from their 
q-bosonic counterparts by replacing q by —q. The unsigned q-Lah numbers 
are met in the normal ordering of ((A^) r A s ) n for r = 2 and s = 1 [15] and 
their q-analogues have been obtained in [4,5]. 

4. Anti-Normal Ordering of q-fermion operators 

In this section, we seek an expansion for (AA^) r . First, we consider q- 
bosonic operators in (3). It is straightforward to expand 

(AA^Y = Y,A:A s {rf)% (22) 

s=l 

and use 

(At)M = lA(At)--i[ a ](At)'-\ 

to obtain a recurrence relation for A r s with A\ = 1, A r = 0, A r r+1 = as 

A r s +1 = q-^AU-isjqSAl. (23) 

A similar relation for q-fermionic operators anti-normal orderering, namely 

(FF^Y = ^^F s (F f ) s , (24) 

s=l 

can be obtained from (7). We have from (7) 

(F^) S F = (-l) s q' s F(F^ s - (-l) a g-'[s] / (F t ) s - 1 . 
Using this and (24), we find (B\ = 1; B r = 0; B r r+1 = 0), 

B r s +1 = (-ly-'q-^BU-i-iyq-^fBl. (25) 
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These recurrence relations (23) and (25) are different from those of q-Stirling 
or Lah numbers. 

In order to obtain a relationship between A r s (B r s ) and ([71]/), we use 
the q-boson Fock space states for (22) and q-fermion Fock space for (24). 
Then using the standard action of the creation and annihilation q-operators 
on the Fock space states, we find 

[B+1B = ±a: 1 ^, 

[«+!]/ = tf-lffi' <26) 

It is interesting to observe that while the normal ordering of operators yielded 
expressions for [n] r b ([n]j) in terms of q- Stirling numbers of the second kind as 
'falling factorials' (namely (14) and its q-bosonic analogue), the anti-normal 
ordering yields expressions for [n]l([n] r f) in terms of A r s , B r s as 'raising fac- 
torials '. 

5.q-fermionic Stirling number of second kind - a probabilistic view 

In this section we extend the theory of product densities of Ramakrishnan 
[10] to a q-extension of the stochastic variable [n(E)]f which depend on a 
continuous parameter E taken to be ordinary variable, The statistical prop- 
erties of these q-fermionic stochastic variable taking values [n(E)]f will be 
governed by q-stochastic point processes. Recall that [n]f for q < 1 never ex- 
ceeds unity. So the feature assumed in [10] namely atmost one particle occurs 
in the interval dE is maintained. Now we define the q-number of particles in 
the range E and E + dE to be [n(E + dE) — n(E)]f which is just [dn(E)]f. 
Following [10], we take the probability that there occurs [1]/ = 1 particle 
in the interval dE is proportional to dE and that for the occurrence of [n]f 
particles in dE is proportional to (dE) n . The average number of particles 
in the interval dE, denoted by £([dn(E)]f), is represented by a q-function 
f[ q \E) such that 

£([dn(E)]f) = f[ q \E)dE. (27) 
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Denoting the probability that [n]f particles occuring in the interval dE by 

P q ([n]f), we have 

P q ([l]f) = PoW = A q \E)dE + 0{{dEf), 

P q ([0]f) = PM = l-f[ q \E)dE-0{{dE)% 

P q ([n] f ) = 0((dE) n ) 5 n > 1. (28) 

The average of the r th moment of [n]f is then 

£(M/) = ENW4 

n 

= f[ q \E)dE = E{[dn{E)\), (29) 

where the second step follows from (27) and (28). Thus all the moments 
are equal to the probability that the q-stochastic variable assumes the value 
[1]/ = 1. This feature of [10] is maintained here. f[ q \E) is the q-product 
density of degree 1. 

Now we consider the distribution of [n]f particles in the E-axis, that is, 
in the intervals dE\, dE 2 , ■ ■ ■ dE n , with [l]/(= 1) particle in each interval. For 
the first interval dE\, this can be done in ^ / C[i] / = — 1]/! = [n]f 

number of ways, where ' n ^C[i] / is the q-binomial coefficient. Thus, from the 
average number of particles in (27), we have 

f[ q \E 1 )dE 1 = [n] f f[ (l)0 (E 1 )dE 1 , (30) 

with 

/ f[ q) \E)dE = 1. (31) 

J whole range 

If we now use the reamaining particles as ([n]f — 1), then the number of ways 
of putting [1]/ particle in dE 2 will be [n]f — 1. Then the joint probability of 
putting [1]/ particle each in dE\ and dE 2 will be proportional to [n]/([n]/ — 1), 
if we were to use the same product density in [10]. But, in this way, we will 
not be exhausting the total number [n]f of particles. So, in dealing with q- 
numbers, it is necessary to introduce q-product densities, such that the joint 
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probability of putting [1]/ particle each in dE\ and dE 2 will be taken to be 
proportional to [71] /[Vi — 1]/- Then the average number in this case will be 

S([dn(E 1 )] f [dn(E 2 )] f ) = f^ q) (E 1 , E 2 )dE 1 dE 2 , 

= [n] f [n - l]//i (9)0 (Si)/i (9)0 (^) dE x dE 2 . (32) 



Proceeding further, the joint probability of putting [1]/ particle each in 
dEi,dE 2 , ■ ■ ■ dE n will be proprtional to [n]f\, thereby exhausting the total 
number [n]f particles. This gives 

f^(E 1 ,...E m )dE 1 ...dE m = W/! & ) \E 1 )...f[ q) \E m )dE 1 ...dE m 

[n - m\ f \ 

tf\E 1 ...E n )dE 1 ...dE n = [n\f\f[ q) \E 1 )..-f[ q) \E n )dE 1 ---dE n . (33) 

In (32) and (33), the intervals do not overlap. When the intervals overlap, a 
degeneracy occurs [10] and then for a finite interval AE = E u — E t 



SddniE.MdniE^) = f[ q) {E)dE 

t J Et J Et 

rE u rE u . . 

+ / / f 2 iq) (E u E 2 )dE 1 dE 2 .(34) 

JE t JE t 



The r th moment of the q-number of particles in the finite range AE, namely 
£([ n ]/ ae) can be represented, after taking the degeneracy into account, by 



*(["]/ a*) = E^; ■■■ r f( q \E 1 ,..-E s )dE 1 ...dE s , (35) 

JEt JEt 



where the coefficients C r s are functions of r and s alone. For [n]f fixed, 
integrating over the whole range and using (31) and (33), we obtain 

which is same as (14) upon identifying C' s with T' s . This derivation of (36) 
gives the role of the q-fermionic Strirling number of the second kind as taking 
into account the degenarcaies in the joint probabilities of the distribution of 
q-fermionic number [n]f as a stochastic variable when the intervals overlap. 
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6.Bargmann Space Representation for q-fermion operators 



A Bargmann space realization for q-bosons has been developed by Bracken, 
MacAnally, Zhang and Gould [16] and that for q-fermions has been developed 
by the author [17]. As the q-fermion coherent states involve quasi-Grassmann 
variable vp, the space consists of monomials of ip. In this section, we briefly 
recall the main results to illustrate two points. First, a naive replacement of 
q by — q is not sufficient and is incorrect. Second, we would like to represent 
(9) and (24) in terms of 'differential operators'. The q-fermion coherent state 
[13] is given by 



\*> = (ef^)"V Ft |0>, 



[2n + 1] 



f- 



2n+l>}, (37) 



and it can be verified F\ip >= >. We will first map a vector in the 
Hilbert space \<p > to a function <f)(ip) by 

</#) = <^ f |0>, 

_l oo i2n 

<2n+l|0>}, (38) 



+ 



\l^n + 1] 



f- 



where we have made use of the anti-commuting property of with F and 
F^ . In [17] it has been shown that (j>(ip) is an entire function. 

Now consider the matrix element < ^\F^\(f> >. Using F\ip >= ip\ip >, 
we have 

<^ t |i ?t |0> = WW), (39) 

and so in the space of(p(ip), F^ is represented by multiplication by ip. Consider 
now the expression ip < ^\F\(f) > which can rewritten as 

i)<^\F\(f)> = <^ t |F t F|0>. (40) 
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Writing F^F as [N]f (which is possible in view of the expression for \ip > in 
terms of expansion of q-fermion Fock space states), we have 

ij) <ifi\F\<j>> = <^\[N] f \(f)>, 

= <^ t | 1 " 1 ( ' 1)jVgJV |0>, (41) 



where in the last step we used (8). Using (38), we have < ip^\q N \<j) > = <p(qip) 
with the prefactor in (38) unaltered and so, 

^<^|F|0> = (42) 

This suggests to introduce q-differentiation as 

= m^t^l. ( 43) 



d q V "' ^(i + g) 

Thus in the space of<f>(if)), F is represented by q-differentiation with respect 
to ip. It can be verified j^ip n = [n\fip n ~ 1 . In this way, we realize that the 
Bargmann spaces for q-boson and q-fermion are very different. The mono- 
mials of quasi-Grassman variables cannot be obtained by naive replacement 
of q by —q. Now using these results, the expression (9) can be written as 

- (44) 

and the expression (24) gives 

= t«?(^)V#W. (45) 

Expressions (44) and (45) reveal the role of q-fermion Stirling number of 
the second kind and 'number £>" appearing in the anti-normal ordering of 
q-fermion operators, in expressing powers of q-differential operators. 

7. Summary 

We have considered the q-fermion numbers introduced in the q-fermion 
oscillator algebra by Parthasarathy and Viswanathan [6] in detail and ob- 
tained q-fermionic Stirling numbers of first and second kind explicitly. The 
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q-fermionic Bell number is obtained by means of q-fermionic Dobinsky for- 
mula, q-fermionic Lah numbers are also considered. These results agree with 
those of Schork [4,5]. The case of 'anti- normal ordering' of q-fermionic anni- 
hilation and creation operators is studied and expansion coefficients A r s for 
q-bosonic operators and T r s for q-fermionic operators are introduced. Recur- 
rence relations for these are derived and these are very different from those 
encountered in q-stirling numbers or q-Lah numbers. In this sense they are 
new. By taking matrix elements of the defining relations for these coefficients 
between q-bosonic and q-fermionic Fock space states, we obtain expressions 
for the powers of q-bosonic and q-fermionic numbers in terms of 'raising fac- 
torials'. These compliment the role of the q-Stirling numbers of the second 
kind as they express powers of q-bosonic and q-fermionic numbers in terms 
of 'falling factorials'. The theory of product densities of Ramakrishnan [10] 
is extended to q-stochastic point process and the necessity of introducing q- 
product densities is emphasized. This leads to the identification of the effect 
of the degeneracy with q-fermionic Stirling numbers of the second kind. We 
have given a Bargmann space representation of q-fermion operators F and 
using q-fermion coherent states. This representation is used to express 
powers of q-differential operators acting on the space of entire functions of 
quasi- Grassmann variable as a series. 
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